Relation between Magnitude Series Correlations and Mult ifr act al Spectrum Width 
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We study the correlation properties of long-range correlated time series, Xi, with tunable corre- 
lation exponent and built-in multifractal properties. For the cases we investigate we find that the 
correlation exponent of the magnitude series, \xi+i — Xi\, is a monotonically increasing function of 
the multifractal spectrum width of the original series. 



PACS numbers: 87.10.+e, 87.80.-(-s, 87.90-(-y 

Many natural systems exhibit complex dynamics char- 
acterized by long-range power-law 2-point correlations 
[1-3] . In some cases the output of such systems — a series 
of a fluctuating variable Xi — may be characterized by a 
single 2-point correlation exponent, while in other cases 
the series may be quantified by a spectrum of exponents 
cafled the multifractal (MF) spectrum [4-6]. The MF 
spectrum quantifies in details the long-range correlation 
properties of a series Xi. Nonetheless, additional way to 
measure the correlation properties of a complex series, 
Xi, is to analyze the long-range power-law correlations 
for the magnitudes of the increments \ Axi \ = \xi^i —Xi\. 

For example, complex fiuctuations in financial indices 
are characterized by a broad MF spectrum [7]. Moreover, 
recent reports showed that there are positive long-range 
power-law correlations in the volatility (magnitude of the 
increments) in such financial time series [8]. Biomedical 
time series are also found to exhibit similar complex be- 
havior — e.g., the human heartbeat dynamics generates 
fiuctuations which are anti-correlated on a wide range 
of time scales [9] and are also characterized by a broad 
MF spectrum [10]. Furthermore, it was recently reported 
that the magnitudes of these fiuctuations are positively 
correlated and exhibit power-law scaling behavior over 
scales from seconds to hours [11]. These and other ex- 
amples [12,13] raise the possibility that there may be a 
relation between the scaling in the magnitude series and 
the MF spectrum of the original series which is indepen- 
dent of the type of long-range correlations present in the 
original series. 

Here we study the relation between the scaling expo- 
nent which quantifies the power-law correlations in the 
magnitude series and the MF spectrum of artificially gen- 
erated time series with analytically built-in MF proper- 
ties. We find that the scaling exponent of the magnitude 
series is a monotonically increasing function of the MF 
spectrum width of the original series. 

Long-range power-law correlated series are monofrac- 
tal when they are characterized by exponents T{q) for 
different moments q being linearly depend on q, 



with a single Hurst exponent [14], 

H = dr/dq = const. 



(2) 



Other long-range correlated series are multifractal and 
are described by T(q) that depend nonlinearly on q, 



h{q) = dr/dq ^ const. 



(3) 



where h{q) is the local Hurst exponent [6]. The MF spec- 
trum of the series is defined by 

D{h) = h{q)q-T{q), (4) 

where for monofractal series it collapses to a single point. 



D{H) = -t{Q). 



(5) 



Here we show that the magnitude series scaling exponent 
is related to the width of the MF spectrum D(h) of the 
original series. 



(6) 



T{q)=Hq + T{ld), 



(1) 



where /i„im and hmax are the minimum and maximum 
of the local Hurst exponent. 

We use the algorithm proposed in [15] to generate 
artificial MF series. Briefiy, the random MF series is 
built by specifying discrete wavelet coefficients which are 
constructed recursively using a random variable W (see 
[15,16] for details). Once the wavelet coefficients are ob- 
tained, we apply inverse wavelet transform to generate 
the MF random series; we use the 10-tap Daubechies 
discrete wavelet transform [16,17]. The MF properties of 
the obtained artificial series depend on the form of dis- 
tribution from which the random variable W is chosen. 

We consider two different types of probability distri- 
butions for the random variable W , the log-normal dis- 
tribution and the log-Poisson distribution. For these two 
examples the MF properties are known analytically [15]. 
The MF spectrum D{h) is symmetric for the log-normal 
distribution and asymmetric for the log-Poisson distribu- 
tion. 
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First wc consider the case of a log-normal random vari- 
able W (Fig. la) — i.e., ln|W^| is normally distributed 
and fj,, are the mean and the variance of ln|iy|. We 
choose n = — (ln2)/4. In this case the scaling exponents, 
r(g), are given by [15] 



f \ <^ 2 , 9 1 



(7) 



and the MF spectrum D{h) is 



Dih) = - 



(/t- 1/4)^ In 2 
2^ 



+ 1. 



(8) 



We calculate the MF spectrum width by solving D{h) 
and find 



Aft 




(9) 



For (7 = 0, the scries is monofractal, while for large value 
of a it has a broad MF spectrum. The MF spectrum is 
symmetric around D(l/4) = 1 and collapses to a single 
point D(l/4) = 1 when cr ^ [Eq. (8) and Fig. 2a,b]. 

We generate series with different a values ranging from 
to 0.1 (20 realizations for each cr) and calculate the mag- 
nitude scaling exponent using the method of detrended 
fluctuation analysis (DFA) [18] (see Fig. Ic and Fig. 3). 
Wc find that the scaling exponent of the magnitude scries 
increases monotonically with a from uncorrelated mag- 
nitude series with scaling exponent w 0.5 for cr = to 
strongly correlated magnitude series with exponent « 1 
for (7 = 0.1. Further, the magnitude exponent converges 
to 1. 

Next we show that the scaling exponent of the magni- 
tude does not only depend on the positive moments (and 
especially the fourth moment) but rather relates to the 
entire MF spectrum. For this purpose we study an ex- 
ample with an asymmetric MF spectrum where the expo- 
nents for negative moments, h{q < 0), are changed dras- 
tically when the MF spectrum width is changed, while 
the exponents for the positive moments h{q > 0) are less 
significantly changed. We tune the parameters in such a 
way that the fourth moment t(4) is fixed and the second 
moment t(2) is hardly changed. 

To generate a series with an asymmetric MF spectrum 
we consider the case of a random variable W from a log- 
Poisson distribution, where A is the mean and the vari- 
ance of a Poisson distributed variable P, and ln|Ty| has 
the same distribution function as P In ^ -|- 7. We choose 
A = ln2 and 7 = — (5'*(ln2)/4 where S is an appropri- 
ately chosen positive parameter [15]. Under this choice 
of parameters. 



and 



D{h) 
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(10) 
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Thus the fourth moment t(4) is independent of 6 

r(4) = 0, (12) 
and the width of the MF spectrum depends on 5 

Ah = ((5V4 -elnS)- (6^/4) = -e\nS. (13) 

For 5=1 the MF spectrum collapses to a single 
point (monofractal) with Hurst exponent H = 1/A and 
D(l/4) = 1. For both cases, decreasing or increasing 5, 
the MF spectrum becomes broader. Here we change the 
MF spectrum width by decreasing 5. The major change 
in the MF spectrum occurs for negative moments q < Q 
(Fig. 2 c, and d). We find that the scaling exponent 
of the magnitude series increases monotonically with the 
MF spectrum width (Fig. 3) although the positive mo- 
ments are hardly changed. Thus, the magnitude series 
scaling exponent is related to the entire MF spectrum of 
the original series and not just to the positive moments. 

We summarize the results of the log-normal and the 
log-Poisson examples in Fig. 3 the magnitude scaling 
exponent is plotted versus the width of the MF spectrum 
[Eqs. (9), (13)]. Surprisingly, these two examples collapse 
on the same curve. This collapse suggests a possible one- 
to-one relation between the magnitude scaling exponent 
and the MF spectrum width [19]. 

The relation between the scaling exponent of the mag- 
nitude series and the MF spectrum width of the original 
series is consistent with recent works on MF random walk 
models with built-in correlations [20]. Recently [11] it 
was found that power-law correlations in the magnitude 
scries are related to nonlinear [21] properties of the in- 
crement series Axj — nonlinear series have correlated 
magnitude series while linear series have uncorrelated 
magnitude series. On the other hand, monofractal series 
are linear series while MF series series are nonlinear [10]. 
Thus the origin of the magnitude series correlations and 
the corresponding MF spectrum may be related through 
the nonlinearity of the original time series. 

The results described here may be important from 
a practical point of view. The calculation of the MF 
spectrum from a time series involves advanced numerical 
techniques [6] and requires long time series. Our analysis 
of the magnitude series is less sophisticated than the MF 
analysis and is applicable to shorter time series than one 
need for the MF analysis. We note that a direct numer- 
ical calculation of the MF spectrum width from a time 
series may not follow the monotonic relation described in 
Fig. 3 due to overestimation of the MF spectrum width 
caused by (i) the numerical technique for calculating the 
MF spectrum [and estimation of T{q ioo)] and by 
(ii) finite series length. However, we expect the magni- 
tude scaling exponent to increase monotonically with the 
width of the MF spectrum [12,19] since both are related 
to the nonlinearity of the original series. 
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FIG. 1. (a) An example of 32,768 data points of an artifi- 
cial MF random series with Hurst exponent oi H = 0.25. We 
use the log-normal W wavelet CEiscade with /j = — (ln2)/4 and 
a = 0.1 to generate the artificial series, (b) The magnitudes 
of increments obtained from the original scries shown in (a). 
Patches of more "volatile" increments with large magnitude 
are followed by patches of less volatile increments with small 
magnitude. In contrast, a monofractal series in (a) would re- 
sult in a homogeneous series for the magnitudes of increments 
|Aa::i[. (c) The root mean square fiuctuations F{n) versus the 
window size n (obtained by the second order detrended fluc- 
tuation analysis [18]) for the original series shown in (a) and 
for the magnitude series shown in (b). The scaling exponents 
are calculated for window scales n > 64. The correlations 
of the original series follow a scaling law F(n) ~ n^'^^ and 
is consistent with a theoretical Hurst exponent oi H = 0.25 
since H = a — 1 (the difference of 1 is due to integration in 
the DFA procedure). The magnitude series is also correlated 
with an exponent of a w 1. 
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FIG. 2. (a) Theoretical scaling exponent r(g) for cr = 
and a = 0.1 for a MF series generated using random variable 
W from a log-normal distribution. The series has symmetric 
MF spectra where the width of the MF spectrum is propor- 
tional to a [Eq. (9)]. For cr = T(g) is linear and the series is 
monofractal while for a = 0.1 T{q) is curved and the series is 
MF. (b) The MF spectrum, D{h), obtained from T{q) shown 
in (a). For cr = (•) the D{h) MF spectrum is a single point 
indicating monofractal behavior with a global Hurst exponent 
H = 0.25. For a = 0.1 (curved solid line) the D{h) spectrum 
is wide indicating MF behavior. The MF spectrum is sym- 
metric and centered around h — 0.25. The vertical dashed 
lines indicate the minimum (hmin) and maximum (hmax) val- 
ues of the local Hurst exponent, (c) Same as (a) for a MF 
series generated using random variable W from a log-Poisson 
distribution. In this case the MF spectrum is asymmetric for 
which the second moment t(2) and the forth moment t(4) 
are almost fixed. The r(q) spectrum is given for 5—1 (linear 
dependence that indicates monofractality) and for 5 — 0.9 
(curved line that indicates multifractality). Note that T{q) 
remains almost unchanged for positive q's while r(q) changes 
significantly for negative q's. (d) The MF spectrum, D{h), for 
the examples shown in (c). The D(h) spectrum for 5 = 1 is 
monofractal (•); the global Hurst exponent is _ff = 0.25. For 
5 = 0.9 (solid curved line) D{h) MF spectrum is broad and 
asymmetric. D{h) is more stretched to the right indicating 
larger changes in the negative moments. 
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FIG. 3. The magnitude scaling exponent versus the an- 
alytical MF spectrum width. Ah = hmax — hmin- For 
each point we generate 20 realizations each of 65,536 points. 
We calculate the magnitude scaling exponent in the range 
64 < n < 4096 using the second order DFA; the average ex- 
ponent is shown (the standard deviation is loss than 0.05). 
Both examples, log-normal distribution (o) and log-Poisson 
distribution (A), exhibit the same behavior of monotonically 
increasing magnitude exponent as a function of the MF spec- 
trum width. We approximate this increase by l/(l-|-e~^^^''). 
In contrast, the value of the 2-point scaling exponent of the 
original series does not depend on the width of the MF spec- 
trum. 
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